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Abstract 

We consider the dependence magnetization vs. field at zero temperature for the spin-i chains in 
which the intersite interactions regularly vary from site to site with a period p. In the limiting case 
when the smallest value of the intersite interactions tends to zero the chain splits into noninteracting 
identical fragments of p sites and the dependence magnetization vs. field can be examined rigorously. 
We demonstrate explicitly the appearance of plateaus in such dependence and discuss a presence of 
the magnetization values m predicted by the condition p ( | — m) = integer [jjj] . We comment on the 
influence of the anisotropy in the interspin interaction on the magnetization profiles. Finally, we show 
how the case of nonzero smallest value of the intersite interactions can be considered. 
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The theoretical study of (quantum) spin chains attracts much attention during last years. On the one 
hand, a number of quasi-one-dimensional magnetic compounds, the properties of which can be reasonably 
described by the one-dimensional quantum spin models, becomes available. On the other hand, quantum 
spin chains should exhibit various interesting properties the examining of which is of great importance from 
the academic point of view. Thus, the analysis of the magnetization processes at low temperatures may 
yield a step-like dependence magnetization vs. field. The latter problem has received a lot of interest in 
a numerous theoretical, numerical and experimental papers concerning a variety of spin chains and ladders 
ill- 

In what follows we discuss one mechanism which generates a step-like dependence of the magnetization 
vs. field, that is a regular alternation of the intersite interactions. Namely, we consider a chain of N — > oo 
spins i governed by the Hamiltonian 

N N 
n— 1 n— 1 

assuming that the intersite (antiferromagnetic) interactions </" (> 0) vary regularly from site to site with a 
period p, i.e. a sequence of parameters is { Jf , Jf, Jf, . . . , J£, Jp,Jp, Jf, Jf , Jf , ■ • •}■ If J^ — J^ = J% = Jn 
Eq. (0) is the Hamiltonian of the isotropic Heisenberg (XXX) chain, if J%,J% ^ 0, J* = Eq. (|]) 
corresponds to the anisotropic XY chain. For the latter chain one can differ two limiting cases, namely, i) 
Jf% = Jn = Jn — the isotropic XY (XX) chain and ii) J* = J n , J% = — the Ising chain. Besides, the 
Hamiltonian (|l]) contains a uniform external field h directed along z axis (that is called the transverse field 
for XY chains). 

We shall be interested in the magnetization per site m = jf^2n=i( s n) (the- angle brackets denote the 
thcrmodynamical canonical average) or more precisely in the dependence m vs. h at zero temperature. 
Oshikawa et al. Q| using the Lieb-Schultz-Mattis theorem and the bosonization techniques for a general 
quantum spin s chain with axial symmetry argue that the magnetization obeys p (s — m) = integer, i.e. in 
the case s — ^ the possible values of magnetization are m = p ~^ fc , k = 0, 1, 2, . . . ,p. A lot of work has been 
done to check the above mentioned criteria using various approximate analytical approaches and numerical 
techniques. On the other hand, for the XX chain in a transverse field the magnetization m can be calculated 
rigorously. The result at zero temperature reads 

1 f°° 

m=-J dEp(E)(29(E) - 1), (2) 

where p(E) is the density of states of the Jordan-Wigner fermions which is known explicitly for any finite 
period p [Q and 0(E) is the Heaviside step function. Regular nonuniformity leads to a splitting of the 
fermion band of the uniform chain into several subbands that, in its turn, immediately leads to plateaus in 



the dependence m vs. h as follows from the above formula (|2|) for m. The values of the characteristic fields 
at which plateaus start and end up are the solutions of two algebraic equations of pth order, whereas the 
possible values of m are connected with the possible differences in the numbers of subbands at E < and at 
E > (see (g)). One more related work concerning the regularly alternating transverse Ising chain has been 
reported in Ref. ||. Contrary to the transverse XX chain, the regular alternation of exchange interactions 
for the transverse Ising chain does not lead to plateaus in the dependence m vs. h. 

In what follows we consider in some detail a particular case of the regularly alternating spin-i chain (|l|) 
when the smallest value of the intersite interactions equals to zero. Without a loss of generality we may put 
Jp — 0. In such limiting case a simple picture for explanation the zero temperature magnetization profiles 
emerges. Really, in this limit the chain consists of the noninteracting clusters every one of which contains 
p sites. The magnetization of the chain per site m follows from the magnetization of the cluster M p after 
dividing by p, and M p = (GS|5 P |GS), where S p = s\ + . . - + s p and |GS) is the ground state eigenvector of the 
cluster Hamiltonian. The appearance of plateaus arises due to the change of the ground state with varying of 
the field. Such a viewpoint is known as the strong-coupling approach. It was exploited in a number of papers 
devoted to the spin chains with a periodic modulation of the intersite interactions and the spin ladders 
We discuss the strong-coupling limit to get a better understanding of the obtained earlier rigorous results for 
the transverse XX and transverse Ising chains by means of the continued fraction approach [Q, || as well as 
to discuss the influence of the anisotropy in spin interaction on the zero temperature magnetization profiles. 
On the other hand, demonstrating how does the strong-coupling approach work in the exactly solvable case 
we can reveal a region of validity of this approximate method. 

We start with the regularly alternating XX chain with p = 2. Assuming J 2 — one splits the chain into 
noninteracting clusters containing two sites. The Hamiltonian of the cluster reads 

H 2 = -h(sl + s z 2 ) + J 1 (sfs% + s y 1 s v 2 ) . (3) 

The eigenvalues of (||) are 

E 1 = ~J 1 , E 2 = -h 7 E 3 = h, Ei = ^J x] (4) 
the corresponding eigenvectors are 

|i> = ^(ITii2>-UiT a », |2) = |TiU |3) = |lx| 2 ), |4> = -±=(|Ti42> + |liT 2 ». (5) 

Moreover, as it follows from (|) (ll^l) = 0, (2|S , 2 |2) = 1, (3|S2|3) = -1, (4|S , 2 |4) =0. For < h< |Ji 
one concludes from (Q) that |GS) = |1) and therefore M 2 = 0. For iJi < h one finds that |GS) = |2) and 
therefore M 2 = 1. Similarly the case h < can be considered. As a result one deduces that the magnetization 



curve to vs. ft should exhibit a plateau at to = (if — |Ji < ft < \J\) and at to = ±i (if ft > |Ji and 
ft < -|Ji). 

Let us consider further the case p — 3. The cluster Hamiltonian is as follows 

h 3 - -ft (4 + 4 + 4) + (44 + 44) + j 2 (44 + 44) . (6) 

The eigenvalues of the Hamiltonian (0) are as follows 

£i = -^ft-i\/^ + Jl E 2 = hi-^Jl + Jl E 3 = ~h, E A = ~h, 

E 5 = —Ei, E 6 = —E 3 , E 7 — —E 2 , Eg = —E x . (7) 

Moreover, (3|5 3 |3) = §, (1|S 3 |1) = <4|S 3 |4) = (7\S 3 \7) = ±, (2|S 3 |2) = (5|S 3 |5) = (8|5 3 |8) = -|, <6|S 3 |6) = 
— |. Therefore, for < ft < | ^/ Jj 2 + Jf (since |GS) = |1) (see (^))) one finds M 3 = | and m = |, whereas 



for \\f J\ + J% < ft (since |GS) = |3)) one finds M 3 = | and m = 4. As a result one concludes that the 



dependence m vs. ft exhibits plateaus at to = ^ (if \\/ 4 + 4 < ft), at to = g (if < ft < \\J 4 + 4)i a t 



to = -\ (if -^V^ + J 2 < ft < 0), and at m = -~ (if ft < -|V^i + J 2 2 )- 

It is interesting to note that the relevant for zero temperature magnetization low-lying levels of the 
Hamiltonians (||) and (||) follow from the following Hamiltonians 



1 



H 2 = -hS z + ^J 1 [(S z f-l), S* = {±1,0} (8) 



and 



W3 = -ft^ + J//^(W-~), 5* = {±|,±i|, (9) 

respectively. The appearance of plateaus becomes evident from ^ (or (^|)) since, e.g., at small ft > spins 
are fixed to S z = (or to S z = \) and for large ft > they should be in the S z = | (or 6> 2 = 1) state. 
We pass to the case p = 4. The eigenvalues of the cluster Hamiltonian are as follows 



E 1 = --\/Ji + (j 1 + j 3 y = -E 167 



1 



E 2 = -h- ^\I4 + 4 + 4 + yJf + 4 + 4 - 2 44 + 2 44 + 2 44 = - E m, 



i 



e 3 = ft- ^\I4 + 4 + 4 + \I4 + 4 + 4 - 144 + ^44 + 244 = -Eu, 



E i = -^Ji-lJ l -J,r- -,E, : . 



1 



E 5 = -h-^= x / j( + Ji + Ji -^4 + 4 + 4 - 2 44 + 2 44 + 244 = - E i2, 



1 



E 6 = h-—^4 + 4 + Ji - ^Jjf + 4 + 4 - 244 + 244 + 244 = -En, 

E 7 = -2h = -E w , E s = = E 9 . (10) 



Usually, a chain with p = 4 (e.g., if J± = J2 = J3 = J) exhibits plateaus at to = for — /ii < h < hi 
when |GS) = |1), (1|5 4 |1) = 0, at m = \ (m = -|) for hi < h < h 2 {-hi < h < -hi) when \GS) = |2), 
(2|S* 4 |2) = 1 (\GS) = |3), (3|5 4 |3> = -1), and at m = \ (m = -|) for h 2 < h {h < -h 2 ) when ICS 1 ) = |7), 
(715*417) = 2 (|G5) = 1 10} , (10| J S , 4 j 10) = —2). However, for special values of parameters Ji, J 2 , J3 not all 
possible values of m, i.e. 0, ±j, ±i, are observed. For example, if Ei = E% in ( |l0| ) (this occurs when 
Ji = J 2 = J, J3 = 0) the plateau at to — disappears. In terms of the Jordan- Wigner fermions in such a 
case one has § p(E) = \5{E + h + ^J) + \5{E + h) + \5{E + h- & J). Therefore, in accordance with (|) 
to = occurs exactly when h = since any small positive (negative) h immediately yields to = ^ (m = — i). 
Let us turn to the transverse Ising chain. For p = 2 the cluster Hamiltonian eigenvalues are 

E 1 = -^h* + j£j? = -Ei, E 2 = -^J 1 = -E 3 . (11) 

Note, that the ground state is |1) for any h. Moreover, ( 1 1 1 1 ) = if h = and (IIS2II) — ► 1 if /i — ► 00. Thus, 
the considered chain does not exhibit plateaus that is in agreement with the result for a general (without 
restriction J 2 = 0) transverse Ising chain with p = 2 reported in ||. In fact, the absence of plateaus in 
magnetization curve is not surprising and is conditioned by the different symmetry of the transverse Ising 
chain. Thus, Y] n s z n for this model with arbitrary J 2 does not commute with the Hamiltonian (in contrast 
to the case of the transverse XX chain) and hence (GS\S P \GS) should vary continuously with changing h. 

One can convinced himself that the formulae (^), (|7|), (|Io|), ( fi"l| ) are valid for ferromagnetic sign of all 
or a part of interactions and hence the described picture is not restricted to the case J„ > 0. For the 
isotropic Heisenberg chain with p = 2 with the antiferromagnetic interaction 3\ > the cluster Hamiltonian 
eigenvalues are Ei — — | Ji, E 2 = —h + iJi, E3 = ~Ji, E4 = h + i Ji, and therefore one concludes that 
m = for —Ji<h<Ji and m = \ (m = — h) for Ji < h (h < — Ji). However, for the ferromagnetic 
interaction Ji < one finds instead Ei = —h— t|Ji|, E 2 = —\\Ji\, -E3 = h— j\Ji\, E4 = — ||Ji|, and hence 
the plateau at to = does not appear. 

Evidently, in the considered limit J p — an arbitrary type of the intersite interaction (the anisotropic 
XY chain, the Heisenberg-Ising (XXZ) chain etc.) and the value of spin s can be examined easily. In 
addition to the total magnetization the on-site magnetization can be calculated in a similar way. Besides, 
the treatment in that limit can be applied to the spin systems of higher dimensions. 

Finally, let us discuss briefly how the obtained results can be used for the analysis beyond the limit 
J p = 0. The transverse XX chain or the transverse Ising chain can be studied rigorously [5) in contrast 
to the Heisenberg or more complicated chains. Different ways to take into account a nonzero value of the 
smallest interaction perturbatively at certain value of m have been elaborated (j| . Let us demonstrate how 
it can be done considering for concreteness the regularly alternating XX chain of period 2 in which now 



Ji > J2 / 0. The Hamiltonian of that system naturally splits into two parts 

H = H Q + V, (12) 



1 



H Q = ...-h (sf 01 + sf 02 ) + - Ji (s^oi 3 ^ + s 



'101 J 102) 



V — ... + -,h (Sl02 s 103 + s 102 s 103) + 

Only the two lowest levels of the two-site cluster discussed above will be taken into account. We assume 
h > 0; the relevant states arc |1) 51 = -±= (| T101I102) - I I101T102)) and |2) 51 = | Tioiti02>- Let us introduce 
new spin \ operators of which act in a following way 

fffl|l>Bl = -^|l>51, <i|l>51 = |2) 5 1, ^|1)51 = 0, 

^12)61 = ^12)51, <7 5 + i|2)5i=0, ^12)51 = 11)51, .... (13) 
Then H can be written approximately (only the lowest levels are reproduced) as 

^0 = t {-\J, Q - of) -h(±+ of)) =-\(h + \ Jr) L-(h- \ Ji) t a t (14) 
with L = ±N. Really 

\(h+ -(h- l -J^j of} |2), = -h\2) u (15) 

that are just the two lowest levels of Ith cluster. To write the intercluster interaction, e.g., V§i^2 = 
\-h ( s io2 s io3 + s i02 s i03) m t erms 01 the operators of let us consider the following equations 

V51,52|l)51 = ^ ^25^3-^= 1 2) 51 = -^-j= JlS^^tx | l) 5 i , 
^51,5212)51 = \jlsXw\ Tl0lll02) = ^28^03-^=11)51 = ^=J2Sj t 03Cr 5 i|2)51 (16) 

(note, that although | T101I102) = ^75 (| 1) 51 + |4) 5 i), we put | T101I102) = ^75 |l)si since we take into account 
only 1 1) 5 i and |2) 5 i). Hence 

^51,52 = 2^= J 2 (sros^Sl + S 103 o '5l) ■ ( 17 ) 



Further, 



V51,52|l)52 = ^/= J 2 (^^f l 3 ^ 52 + CT 5 1 ^7 = ^l 2 ^ 52 )) = ~ \ J2<J 51°52 I !>52 , 



As a result one concludes that 



V51,52|2) 5 2 = ~^~^ j2(J 5l\ ll03Tl04) = -^^20- 5 f iCT52|2)52. (18) 



^51,52 = -\j2 (^52 + (T 51 cr+ ) (19) 



and therefore 

^=-^(^ + i+«- (20) 
1=1 

The Hamiltonian (|l2|), (fbi]), (|2Cj ) describes the uniform spin-i transverse XX chain of L spins 



l (h+ l -J 1 )L-(h- l -J x ) £ a\ - \j 2 £ (ofof +1 + afaf +1 ) . (21) 



i=i i=i 



Finally, acting like while deriving ([l7|), ([H]) one finds the relation between the operators s" and a" 

+ 1+ _ _ 1 _ z _l/l , \ 

■ ■ ■ i s ioi — ~~ °5H s ioi ~ ~"^7f CTsl ' Sl01 — 2 I 2 1751 / ' 

+ _ 1 + - * - 1 * ^ 

■ s io2 - "yf " 51 ' Sl02 ~~ T/f " 51 ' Sl02 — 2 \2 1751 / ' 

Consider now the magnetization which owing to (22) can be written as m — | (h + Y^2i=i( a i)^j ■ Using 
the well-known results for the uniform spin-i transverse XX chain (see, e.g. Q) one concludes that m = 
at < ft < i {J 1 — J 2) and m = A at \ (J\ + J2) < ft. If ft increases from ^ (Ji — J2) to i (Ji + J2) the 



magnetization increases from to \ as ^ — jp arcsin W 1 — 



(2/1- Ji) 



In Fig. 1 we plotted the zero temperature magnetization profiles for the spin-^ transverse XX chain of 
period 2 with J\ = 1 + S, J2 = 1 — 5 for 5 = 1, 0.8, 0.6, 0.4, 0.2, as they follow from 1) the exact formula 
(|2|) H (solid curves) and 2) the approximate Hamiltonian (|2l]) (dotted curves). These plots demonstrate 
how does the strong-coupling approach work as S deviates from 1. 

To summarize, we have reconsidered the zero temperature magnetization processes in the regularly al- 
ternating spin-i XY chains within the frames of the strong-coupling approach. Besides discussing the 
magnetization plateaus in terms of the spins rather than in terms of the Jordan- Wigner fermions we have 
demonstrated to what extent the strong-coupling approximation can reproduce the exact magnetization 
profiles. 

The author is grateful to J. Richter, N. B. Ivanov, T. Krokhmalskii, O. Zaburannyi and V. Derzhko for 
discussions. He thanks J. Richter for hospitality in the Magdeburg University in summer of 2000 when the 
paper was completed. 
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Figure caption 



Fig. 1. rn vs. h for the spin-^ transverse XX chain of period 2 with J\ = 1 + 5, J2 = 1 — 6, 6 = 1 (a), 
6 = 0.8 (b), 5 = 0.6 (c), 8 = 0.4 (d), S — 0.2 (e), 5 = (f) at zero temperature. Solid curves correspond 
to the exact results following from (|^) Q , dotted curves correspond to the results obtained with the help of 
the approximate Hamiltonian (plf). 
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